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POSER ac t 


The Jolly-Seber method of analysing mark-recapture data 
obtained from a simple homogeneous type of population is 
extended to allow for some peculiarities to be present in 
the population. The type of population considered is 
non-homogeneous, consisting of only two age classes, young 
and adult. The young require a fixed number of years, r, in 
which to mature to adult status. It is assumed that only the 
young of the year and the mature adults are catchable, the 
subadults, for one reason or another, are not. This analysis 
assigns to the young a probability that they survive the r 
years to maturity and a different probability that the 
adults survive one year intervals. Maximum likelihood 
estimates for population size parameters and the survival 
probabilities are provided. 

The estimation equations are applied to a set of 
capture-recapture data obtained from a colony of Pacific 
brant (Branta bernicla nigricans) breeding in the Anderson 
River Delta in the Canadian arctic. These data were 
Bowmiectedricomensesino ASisrbysDr<chiewe Barryxot tne 
Canadian Wildlife Service, with the assistance of the 


author. 
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Introduct ion 
For the past several years concern has been mounting for a 
population of Pacific brant (Branta bernicla nigricans) 
(Pacific Waterfowl Flyway Council, 1978) wintering along the 
west coast of this continent. During the past twenty years 
the wintering flocks of brant have dwindled in their 
traditional areas along the coast from British Columbia to 
California. Government agencies managing the Pacific Flyway 
have noted a synchronous increase in brant numbers along the 
Baja coast of Mexico (Pacific Waterfowl Flyway Council, 
1978). With this growing concentration of geese in one area 
the potential for greater exploitation by hunters becomes a 
real and serious problem. It is also not clear whether this 
shift is caused by brant abandoning their northern haunts in 
favour of southern areas or whether the overall population 
is experiencing actual numerical changes among its many 
breeding subpopulations. Another matter of concern is the 
continuing human encroachment into the coastal habitat of 
this species, a pattern which may simply force the geese to 
move out or else threatens their primary food source, the 
beds of brackish water eelgrass (Zostera marina). The 
transportation of heavy oi] along the west coast and the 
inevitable future production of arctic oi] and gas will most 
likely have a drastic effect on the wintering birds and 
their habitat and on the breeding, molting, and staging 
areas Oncmuneanrctte: coast, Towallevidte these concerns a 


management plan for the Pacific Brant has been formulated 
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(Pacific Waterfowl Flyway Council, 1978), setting forth 
objectives for increasing the numbers and habitat and 
adjusting the distribution of this goose. In meeting these 
objectives there is one major obstacle to overcome: 

The management of a species rests 

upon the data base, and it, unfortunately, 

contains many serious gaps. The 

more serious deficiencies concern 

the amount of habitat and how it 

VWs’ Ghanging), and the-distr 1bution, 

survival and recruitment of brant 

on a subpopulation-by-subpopulation 

basis. These gaps must be closed 

rapidly to provide for well founded 

management. 

(Pacific Waterfowl Flyway Council, 1978) 
hean-year banding *prodr am “ofea’ subpopulation tof *Pacirtic 
brant conducted by Dr. 7. W. Barry with the assistance of 
the author, provides the data base for closing some of the 
Survival and recruitment ‘gaps’ referred to above. However, 
due to a peculiarity in the life history of brant, and geese 
in general, existing methods for estimating population 
parameters are unable to handle this data set. Briefly, the 
pectiittarity TS5eas*tol lows seb rant! general iy do eno tepreed 
Ati d weneymané sinstheip thirdtseason parry? 196/). The one- 
and two-year-old birds do not usually spend the summer with 


the breeding adults, and are therefore absent from the 
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population during banding operations. The goal of this 
thesis will be to develop an analysis technique applicable 
to this particular situation by extending the methods and 
ideas of other researchers in this field. I will then use 
these procedures to investigate the specific subpopulation 
of brant referred to above. 

IT will now discuss the historical eve eomeut OF 
methods for analyzing banding returns, and capture-recapture 
analysis in general, to obtain estimates of various 
population parameters. 

Most developments in capture-recapture analysis 
occurred over the past four decades. The methods dating from 
1930 to 1950, and even into the next’ two decades, are quite 
simple in concept and were designed primarily to estimate 
the size of a population, a useful statistic with which to 
document the status changes of a population. Recently, the 
scope of these analysis procedures has expanded to include 
the estimation of survival and recruitment rates, parameters 
which further the resource manager’s understanding of 
PeSuUhatlon Tuc tuations. 

For the last forty years two major approaches have been 
explored, with varying success, for describing population 
eharacteristics, the deterministic; and. the probabilistic 
[stochastic jlscahalvsiseof <capturesrecaptune data useber., 

1973; Schemnitz, 1980). the first method-involves 
Construct ingedsdeterministic lifeslablesas isadone.im human 


demography. A cohort of individuals of the same age, usually 
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young, is selected from the population and the number of 
these surviving to each age is recorded as this group 
‘moves’ through the different age classes. With these data 
Ginem population *characteni sti.cs@fom eacnrage class. Such as 
the mortality rate or the mean expectation of life, may be 
estimated. This deterministic method, however, is 

Untea istie, 1 gnoming the stochastic Filictuationsrin 
populations and the variability of samples. 

Stochastic analysis, on the other hand, deals directly 
with these features by assigning probabilities to the 
movements between the age classes. This approach has proven. 
to be extremly flexible and has undergone a great deal of 
evolution, resulting in some quite useful, although very 
data-hungry, techniques. 

Stochastic methods of parameter estimation are based on 
information collected by means of some variation or 
refinement of the basic capture-recapture method employed by 
Lincoln in 1930 to estimate the total number of North 
American ducks, N (Seber, 1973). His procedure involved 
selecting a sample of size n from the population, marking 
them in a suitable fashion, and then releasing them back 
into the population. With presumably sufficient time allowed 
for these marked birds to thoroughly mix with the unmarked, 
a second sample of size s is drawn by hunters and the number 
Gesmanked individuals om insthat collectton*noted:. He 
reasoned. thaw the: propor? i on™“or "the (banded birds that were 


shot, and their bands returned, would be approximately the 
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same as the proportion of all of the harvested ducks taken 
frometne tote copulation Thus Be le aid (m/n) = 
(s/N). This equality provides the simple population estimate 
N= 5 x (n/m), referred to as the Peterson Estimate or 
Lincoln Index, the form of which appears throughout much of 
the later literature (Seber, 1973). 

A simple extension of Lincoln’s capture-recapture 
scheme was devised by Schnabel in 1938 (Seber, 1973). He 
allows the sampling to be repeated periodically (for 
example, yearly) for an indefinite time and requires that in 
each sample any unmarked individuals captured be tagged and 
the number recaptured, of those previously marked, noted. 
All of the sampled anew are then released. With both 
Lincoln’s and Schnabel’s sampling plans an estimate of the 
population size only is available. Furthermore, only if the 
population is closed, that is, if death and immigration do 
not occur, or, if the time between samples is instantaneous 
with respect to these features, will these estimates of 
Mimecolm and schmabe! be Valid. The occurrence, of -déatnmand 
immigration alters the proportion of marked {ie WARES in 
Thies population, causing am apparent decrease or an apparent 
increase in the estimate of the population size, depending 
on whether there is a respective increase or decrease in the 
marked proportion. Obviously, this closure requirement 
severely limits the usefulness of early attempts at 
stochastic estimation for the same reason that life tabie 


methods did not work: population levels tend to fluctuate 
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over time. 

Darroch (1959) first broke through the limitations of 
assuming a closed population when he developed procedures to 
account for the possibility of either death or immigration, 
though not both. He refined Schnabel’s long-term census by 
requiring that each captured individual be marked with a 
unique number. Since the complete capture history of any 
individual could now be determined, and hence so could the 
proportion of individuals which were last captured one, two, 
three and so on, years ago be determined, the probability of 
an individual surviving from one sampling period to the next 
could be introduced to the analysis. Darroch showed that, 
with the information retrieved from the recaptured animals, 
estimates of such population parameters as size, survival 
probability, and the probability of capture can be derived 
by maximum likelihood techniques. 

In 1962 Seber extended Darroch’s earlier work by 
estimating parameters of a population affected by both death 
and immigration at the same time (Seber, 1962). One major 
difference, however, existed between Darroch’s and Seber’s 
methods. Darroch collected his information from recaptured 
individuals that were returned to the population alive and 
eouid thus be recaptured at any time sin the future. Seber, 
on the other hand, used data obtained from animals recovered 
by some process which removed them permanently from the 
population fenc. NUN ing =, | ChUs. limi Gindetnes amount .oT. 


potential information available from each marked individual 
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(Seber, 1962; 1970; 1971). Seber’s analysis, however, laid 
the ground-work for further development of a special case of 
parameter estimation by capture-recapture analysis, that in 
which the individuals are recovered only once. A large body 
of literature has accumulated for methods of this type, 
which are extremely useful in the study of heavily hunted 
species of waterfowl and fish. The most comprehensive and 
complete summary is the work by Brownie et al. (1978), which 
describes in detail fourteen different procedures for 
analyzing recovery information based on a variety of 
assumptions. Because these methods were designed for hunted 
species, they are not applicable to the more endangered or 
bess uted ized animals. 

Jolly and Seber, working independently, applied 
Darroch’s now standard Schnabel type census to an open 
population, deriving a general probability distribution 
designed to work with most single populations which are not 
split into isolated sub-populations (Jolly, 1965; Seber, 
1965). They were the first to organize analysis as a concise 
method of achieving max imum likelihood estimates so as to 
account, fors the’ possibility of death, recruitment: and 
immigration of new arrivals into the population, and 
permanent emigration all occuring simultaneously. If the 
sample is taken from a homogeneous population, one in which 
all individuals behave in a similar manner, or equivalently, 
all have the same probability of surviving, then. ,the maximum 


likelihood estimates of survival probabilities, population 
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size, number of recruits, et cetera, have very simple forms 
with obvious intuitive interpretations. Known as the 
Jolly-Seber method, this method of capture-recapture 
analysis was the basis of the more sophisticated techniques 
of estimation to be developed later. 

Because most species do not form a homogeneous 
population and are, in fact, composed of any different age 
groups, for example, young, subadult, and adult, it was 
necessary to make the Jolly-Seber method more realistic and 
flexible by allowing for different survival rates among the 
Various age classes, in other words to create a stochastic 
analogue approximating the rigid deterministic techniques 
mentioned earlier. To this end Stokes (1980) relaxed the 
Jolly-Seber assumption of an homogeneous population while 
using their method to investigate a population of American 
Woodcocks (Philohela minor), a species consisting of two 
recognizable age classes, young and adult. The young birds 
of this species mature to breeding adults in one year, a 
fact Stokes accounted for by assigning different survival 
probabilities to the two age groups. Her methods were 
recently generalized to a population consisting of an 
arbitrary number of one year age classes (Pollock, in 
press). As with all of these models, Pollock's 
generalization provides intuitively appealing maximum 
likelihood estimates for the population parameters. 

This then describes advancements in capture-recapture 


analysis to the point where my analysis begins. The 
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difference between Pollock’s general analysis and my problem 
VSethis: here tf ceai with <a situations ime which the ‘conort of 
young do not return to the nesting grounds until they are 
adult breeding birds. The situation is unique and therefore 


presents special problems to the statistician. 


Xv 


AckKnow ledgments 

I would like to express my deep gratitude to my parents 
who introduced me to the Anderson River country and who 
helped me discover its geese and see its many other charms 
and subtleties. Valuable helicopter and logistic support was 
provided by the Polar Continental Shelf Project of the 
Department of Energy, Mines and Resources. Dr. R. Rout ledge 
of Simon Fraser University found ways to resolve some of the 
difficult obstacles encountered in the derivation of my 
estimation equations. Above all 1 wish to thank Catherine 
McCune-Barry for her patience and her gratefully recieved 


msusoer during this long sproject. 


XV 


wet 


Pies: Sahhii’ et ce 
> no | danke pi muoons itn aa 
snare? frst so. gels dyad: .2hehseype inal bem mo 

beveivet yiluioleio serving: Soria! ton sad OF al 
et ore! eft pride — 


<i 
te 


Table of Contents 


Chapter Page 
LRN eve ap G12 =) mene at ec le So NN xV47I 

Ee Sud pM OUD Sr. nyc. 2° oy eerie mmr een a Cwety ec XV ida 

Ne aE OC Sr Re memr tare Mee eat aris chs 0.06, laliaaen ck a MURR AoA an a ae eae 1 

PRU SSC MD OMmOt oO MOOS lec.) cate eae pre eso eae 1 

ASS UMD Ol) S anit as Geo) 3 tas earn SOR eRe Oa eye ce anne cre 3 

Ney a-gas Nola Race seen ar RMRerry Aly Swart, Siete ae ne eam te S 

Bewleve lopment «Ore NOCH a: tee eee ee Dameron tt sa oe ace amas i, 
ConditiconaleProbabi lity istributions.- +o. 0s i 


VOID ee rOoDado ity DiS Chou L 1 OmerunGc com weet © 


CPSC Ae i cna s pp ee NERS Rel glow Mame men Shek itn ets oe Ree is) 

LU EU arse We HY Zoned oc Mine ee eRe gE Fm yk eer see eR 28 

ee etc AGHRO LINC gta ee ese amyl 1 een eee Ron gies ee ee eee 28 

Bee TOG San mre rer Leta ncomteecc. cy ree ronal eee aie e eae ene ane 30 

CoRR SUPE S een ter arta Mact oer ratte ae terat otic adet an on Mabe: Semmens 32 

Die ONGC US'S Of en eer age ae cma ren) 0 haar Sen ne Pa Retire ees eee Sie) 
(Gelalet O elg tae ee eee On a eek eRe Nr mee as ea ay ee 64 
Cermak Up emiG isl. i025 sie er vk cnents Gc aetna both Soma ieee 70 


Appendix: Derivation of Estimation 
Holey. Me pee ee ee ee tary or ma eR ho we ey Rte i3 


XVI 


v. 


int 
| 2 
_ . 


t- 


7 
7 


> 7) : 
an aan 

: ) 

2 


i. ree: Ss ae Z| - +o tard “7 


eRe pF? Xe ert ; 


des Ee, 
E ae as a ee 5S “att +5 Payee a. 
a. AUPE eee ce SW ys. 2fof tamdaph 
Be MAS. OS Boos peeks a. HORAN 
e8., 8 2) ed Pe [igbot +4 Josmaofeved .8 
enol sugt<terd Vinh Paedeas Lene? Brbed 


a] , : 
el se 


&t... qerfodus rar ut itpne AFT teisto7° are 

rca dykeir’ « «ia oh oneIeues aorlsmetes .2 * 
ee RS ees ee ee 5 6 oe ziayl ead. sis fi 5 
| TCT pnuepisEe A 
OF ewe ape eee she, chsh ois... Sponge 2 
Oe Re er nd 255 ee ee 1, 30 ae eee 


= ae notdsusetT .g 


Cnt fe lehndt 24; Me ig St ers , Re) smitensdi 


rietisnit2? tm Cet se) 60 ; kr erieids 
et) ....- ha teh lied Ce artot Saupe 


List of Tables 


Table Description Page 
I Bancina liSstimibution oT Siti. .j) ST 
12 BandinaeUrstmibuuiOM Of S\ 1 ajyy) 40 
ed Pe CanainGe. 1 Stra OUT On Ors (45 1 2mesy a 43 
IV Summary of Banding Results 46 
V Parameter Estimates 49 
VI Clutch Sizes, Anderson River Delta D2 
Vil Distribution of Replaced Bands oe) 
Vito DiIstmibublioneor Sanding Locations, 

Anderson River Delta 54 
IX Use of Molting Sites by Brant, 

Anderson River Delta 56 


XxV11 


» 
oy 
» 
a 
2 
he 


bt me* EU 


‘erMtibess| fc not uel aheed 
2) (pS1&C yevth sae7Sib 
dnb wBoye aett2 onl 2 fat Fo: eal, 

s7'90 7s0°F AGS yEDAs 


List of Figures 


Figure Page 
die Map Showing the Location of 

the Anderson River Delta 26 
figs - Map of the Anderson River Delta 20, 


xvili 


7 r= _ 


| ss 
D>, See leer Fr 
sited Bgvit rez vabint) si io ant. 


' » 


I. Model 


A. Description of Model 

As already mentioned, Pollock's method for parameter 
estimation in a general age group population is not 
applicable to the Pacific brant population which I will be 
examining. His method allows each one year age group in the 
population to have their own yearly survival rate while the 
Pacific brant require one survival rate to cover the three 
years prior to maturity, and a yearly rate thereafter. To 
amend this deficiency I will derive estimates for a general 
population which is composed of only two discernible age 
classes: young and adult. Members of this population follow 
the simple life history pattern in which the young become 


/ 


full breeding adults after a period of ‘r’ years or seasons 
As@cubadu hess and then maintatmethat adult estatus sunt i | 
death#eror -therdirativon of this maturing phase ‘during each 
breeding season, most of these non-breeding subadults will 
congregate in areas physically separated from the breeding 
portion of the population.’ those “that donot "seperate are 
not easily captured during sampling because non-breeding 
brant have finished molting and are flying before banding 
operations take place. Most mature adults return annually to 
the breeding grounds. The following mode! is based on the 
proceouresrore data collection discussed by Seber (1973) in 


which samples of the population are taken periodically. Each 


unmarked individual captured is marked with a uniquely 
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numbered. band. The number of each recaptured band is 
recorded. At the time of any sample the population is 
divided among four classes of individuals: young, unmarked 
adults, adults which were originally banded as adults, and 
finally, those adults originally banded as young. As the 
banded young mature they enter the latter class of adults, 
while the unbanded young and the adults recruited from 
outside the population feed the group of unmarked adults. 
Further, as the unmarked adults are captured and banded they 
move into the group of adults originally marked as adults. 
This describes the general dynamics between the unmarked and 
marked portions of the population. With this type of 
sampling the capture history, or the times of recapture, of 
any individual is established completely, a feature allowing 
the much finer movements of banded individuals to be traced. 
The treatment of these finer movements is as follows. The 
two types of marked adults, those banded-as-adults and those 
banded-as-young, are further grouped according to their 
capture histories, one group corresponding to each possible 
history. The unbanded adults form the group with no capture 
history. At each sampling stage, depending on whether they 
are captured or not, the members of each capture-history 
group can be shunted from one group to another. This 
ateatysise tUunetioms by documenting these exchanges 
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Assumpt ions 

The following assumptions, typical of those originally 
made by Jolly and Seber and used in later extensions of 
their method, appear to be necessary for this type of 
analysis to be tractable. I have modified and added to them 
SOmunNclatNCy altatmewDan ticular $1 (uation id tl eam 
considering. 

1. Every animal in the population, whether marked or 
unmarked, behaves independently of the others and has the 
same probability of being captured in the i’ th sampie. This 
probability is conditional on the individual being alive and 
present in the population when the (i+1)'th sample is taken. 

2. Every adult has the same probability of surviving 
from the i’ th to the (i+1)’ th sample and of being in the 
population at the time of the (i+1}’ th sample, given that it 
was alive and present in the population immediately after 
the i’ th sample is released. 

3. Every young individual has the same probability of 
surviving from the i’ th release to the (itr)’ th sample and 
of being in the population at the time of the (itr)’th 
sample given that it was alive immediately after the i’ th 
sample. ‘r’ is defined below. | 

4, No deaths occur as a result of the sampling process. 

5. Marked animals do not lose their tags. 

Notation 
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its symbol. 


Time at which the i’ th sample is taken for 


I-40 it cael ewnere. IL is the total number of 


samples taken. i is used to identify each 


sample. 


The total number of individuals in the 


popu la lioncdatecame-a. 


The number of young in the population at time i. 


The number of unmarked adults in the 
HOpu lationea tata Men. 
The number of individuals (adults), first marked 


as adults, that are in the population at time i, 


who were last captured at time j<i, and have a 
capture history w (iw is defined below). 

The number of individuals (adults), first marked 
as young, that are in the population at time i, 
who were last captured at time j<i, and have a 
capture history w. 

The total number of individuals captured in the 
i’ th sample. 

The number of young captured in the i’ th sample. 
The number of unmarked adults captured in the 

1° th sample. 
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Those of the N(i,j,y,w) that are captured in the 
i’ th sample. 

This is the capture history vector of length 
bi, sassociated with the Nid. iew) Niis 7, vaws 
Seiwa alice Ar lly Va Wl 

‘w' will consist of (i-1) elements, where the 
Kon entry win l be a "1" “1fo7P tme individual ts 
captured at time K or a “OQ “1f mot captured at 
time K. A general capture history w, associated 
with N(i,j,w) say, will then look like: 
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The specific capture histories used in this analysis are: 
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(w;0), the capture history of length i with w as 


the first (i-1) elements and 0 as the last. 
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( t:* ), where * is a string of zeros of 


appropiate length. This is notation and means 
that at some previous time an individual had a 
capture history of t, and has not been captured 
since. 

The number of recruits joining the population as 
unmarked adults between the i’ th and the 


(j+1)’ th sample. 
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(ne probability that any individual 15 captured 
im the 3 th samole. 

The length of time required for the 

young of the population to mature and become 
DHeESaInosaaulits. Ff 1S ahi xedea mtn, Value stor 
the duration of the sampling program. 

The probability that a young survives from time 
leet Ouetime: (ier), 

The probability that an adult survives from time 
ing OLA Mela ais) 

The number of marked individuals in the 
peopullaiden at time 1. 

The number of marked individuals captured in the 


i’ th sample. 
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Liaise the number of adults who were captured 
in some sample prior to the i’ th one, not 

Cap tuned i mjpune; ji sthyisamp le, but captured 
subsequently. This expression excludes those 


individuals which were subadults at time i, and 


thus were not available for capture at time i. 
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R(i) is the number of individuals released from 
the i’ th sample which are subsequent ly 

recaptured. Excluded are those which were banded 
as young in the i’ th sample. As with Z(i}, these 
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B. Development of Model 


Conditional Probability Distributions 


The following obvious relations describe the division 


of the population into the four classes, young, unmarked 


adults, and the two groups of marked adults, according to 
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Recalling that a much finer partitioning of the 
population exists, due to the availability of the capture 
history, I will now symbolically describe the probabilistic 
movements between these partitions. My goal is to formulate 
a. joint probability distribution of the numbers of | 
individuals in each of the mark-status classes. To do this I 
Will combine conditional multinomial distributions, using 
Eneamuli tpl cation rule: to form the joint prebabi Jity 
distribution. To save space and to avoid some notational 
drmiacuitves, in othe discussion of theyorobabi lity 
CHSeROU tl ONS mr esna lewrite Pixs) wack not 
distinguishing between a random variable and one of its 
possible values. 

Gonsider Tirst the conditional distribution of the 
number of young caught in the i’ th sample, S(i,y), given the 
number of young in the population, Nli,y), and the 
probability of being captured, P/i). There are two possible 
fates confronting a member of N(i,y). He may escape capture 
at time i or be captured, becoming a member of the Sli,y). 


Following the assumptions, the distribution of S(li,y) is: 
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secondly, the distribution of the unmarked adults in 
ther pOpUlata On dtetime ae 1), Nit+i-ou) and anuthe. 1 th 
sample, S(i,o), is also conditional on several quantities. 
These are the number of unmarked adults present at time i, 
N(i,0), the number of unmarked recruits entering the 
population between time i and (i+1), Bli), the probability 
of capture, P(i), and the probability of surviving to time 
tatiie Ali}. The*distribution is ‘found, as before, by 
documenting the possible fates facing the Nli,o). 
A member of N(i,o) may have one of three things occur 
to it between time i and (i+1). 
1) He can be captured in the i’ th sample and 
hence become a member of the S(i,o) and be 
marked. After release he will next appear, if at 
all, as a member of N(kK,i,w) for some k>i. 
2) He may not be captured and thereby become a 
member of N(i+1t,o}. This class has, however been 
diluted by the B‘i) unmarked adults joining the 
population as recruits between time i and (i+). 
3) He may not be captured and then leave the 
population by dying or permanently emigrating. 
These three mutually exclusive possibilities give rise 


to the following multinomial distribution expression: 
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Similarly between times i and (i+1), members of 


N(i,j,w). have three mutually exclusive fates to choose 


from. 
1) Being captured and becoming a member of S(i,j,w). After 
release he will next appear, if at all, as a member of 

~ N(K,i,w) for some k>i. 
2} Not being captured and surviving to time (i+1), that is 


becoming a member of N(it+1,j,v). 
3) Not being captured and then dying or permanent ly 
emigrating from the population. 
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Mesimiplan distribution 15° formed for =the .S\7,3,V,we and 


N(it1,j,y,v), the adults first marked as young: 
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Joint Probability Distribution Function 

Ne wjomml orobaat laity distribution for ali ~of the 
possible options open to the members of the population, 
taken over all of the samples, is formed by the product of 
the ‘conditional multinomial distributions (2) through (7) 
above. This product is taken over all possible capture 
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LOS ctlet eeelsingerelavions (1) ands Zim torsamo lity 


expressions, the joint distribution becomes: 
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C. Estimation Equations 

To justify the estimation procedures used in this 
sect von wil Mbriewryravscuss valiyves 1965)" aoprodach = te 
this problem, an approach which I have adapted to this 
thesis. Strictly speaking, the population quantities 
aioe Nie ay) and Nit+1,4,),Vi arpeunandom variables. 
affected by the random numbers of individuals surviving over 
time intervals prior to time (i+1), and are not fixed 
characteristics (parameters) of the population. As such the 
probability ’diseributrone( 6) is "not “a ‘true whkelihood’ since 
MO tse not Just “[he*probabl ity oT Cblaining ine observations 
eee SW ole we and oO Py ees DUT "Ot Oo tain iag 
both these observations and the N(i+1t,o), N(it+t1,j,v) and 
Nie iey.yv,v) given the parameters Nit), BVi) = Pt) ArT) and 
Y(i). Jolly (1965) reasons that (8) can reasonably be 
regarded as a likelihood function if population variables at 
time i are being estimated since in this case the Nli,o), 
Nie jaw) and Nii, d,y,W) will be parameters of a distribution 
conditional on events up to time i, and so do not have a 
random component. However, when variables prior to time 1 
qrewpeimo estimated, Jolly (91965) Claims that Intuitively etc 
treat these variables as parameters throughout will lead to 
estimates that are close to the true maximum 1ikelihood 
estimates. Following volly s justification, 1 will proceed 
by considering 96) as a true Iikelihood function b, cr a set 
of population parameters. When deriving the estimate of a 


parameter, say U, it is easier to follow the subsequent 
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derivations if the likelihood function (8) is written as the 
product of two functions one containing all of those terms 
Of (3) Which dresincependent of U {call it K) and the other 
containing all those terms dependent on U, F(U) say. 
Thus es eek or (Ue hewmoti Vationitormthis wa isk become 
clear shortly. 

The maximum 1ikelihood estimates of the three 
probabi litres, Pik esi kK) and ¥(K) are*foundtby 
differentiating the natural logarithm of (8) with respect to 
each of them and equating the result to zero. Since the 
probability distribution function is concave upwards (and 
hence its Jogarithm is concave downwards), the solution of 
this equation will locate the maximum of each parameter. 
Thus comsidering (8) as a function of Pik) (the probabi lity 


of capture): 
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Setting this derivative to zero results in the following 


estimate ot the probability of capture: 
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Similarly the logarithm of the likelihood of the 


probability that the young survive to time (ktr) is: 
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from which the survival rate is estimated by: 
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Intuitively, (10) uses as its estimate the proportion of 
those animals initially marked as young at time kK who have 
managed to survive to become adults at time (kK+r}. 

For the adult survival rate, A(kK), the likelihood 
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Set this to zero-and solve for the estimate of A(kK). 
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Intuitively the denominator of (11) is the number of adults, 
whether marked or unmarked, in the population at time k. The 
numerator is the number of adults at time (k+1), after 
removing the young of the year and any adults which could 
not possibly be present at time k but were at time (kK+1). 


The quotient then is the proportion of adults which have 
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survived to time (kK+1)}). 

The remaining maximum ]ikelihood estimates (for 
integer-valued quantities) are found py equating the first 
backwards difference of the logarithm of the likelihood 


function to zero. That is, for some parameter G, 


if DUG) = Kx kG) 

then Tae) se mie) ee Alane 0G) ) 

anomeee Van Gee th eG ie Int liao 
Sears Ga Neen gs)., 


which is independent of K. 
The maximum likelihood estimate then of G is the greatest 
integer less than or equal to ‘g’ such that {\Intl(g)) = 0), 
fiat esaeNiGr ibe iia Goss). FOke this app ication athougn 
I will consider the maximum likelihood estimate to be the 
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From the above discussion the number of recruits B(k) 


willobe estimated as follows. From (8): 
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Since the relative error in equating Nik+i,o) - B(k) + 1 to 
Wek inoe Bikyetore|argeaNtktive) ts) neqligibile. sinc 
estimate for B(k) is the estimated number of unmarked adults 
at time (K+1), with all of those unmarked adults which had 
survived from time kK removed. This leaves an estimate of 
only those adults recruited into the population between 
times kK and (K+1). The more complex derivations that remain 
constitute the Appendix. . 

To round out the discussion of the model, I will give 
the final forms of the estimation equations and provide an 


imuuitive idiscussion ‘of each. 
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First, the number of marked adults present in the 


population at some time k is estimated by the following 


equation: 


eee M,. = 2 (Sue + Sm) + Sm 
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An intuitive explanation of this equation is as follows. To 
estimate the number of marked individuals in the population 
at time kK, M(k), it is necessary to augment the number of 
these actually seen, Sm(K), by an estimate of the number 
that were present but not seen. There were Z(kK) marked 
animals present in the population at the time, but not seen 
then. However, they are seen at some time in the future. 
Hats cuantity is increased by multipiying it by the ratio of 
the number of marked animals Known to exist just after the 
k’ th sampling (S(kK,o) + Sm(k)) to the number of these that 
are subsequently recaptured, R(kK). If the proportion of 
individuals subsequently recaptured is roughly the same for 
both classes, then Z(k)[S(k,o) a smi] Rix estimates the 
number of marked individuals present at time kK but not 
captured then. By adding this term to Sm(kK), one obtains the 
maximum likelihood estimate of M(kK). This equation for M(k) 
is the key to the whole estimation problem since it is 


involved in most of the remaining estimation equations. 


Ath col 


ot ,awollot 26 @F nett aype BetltF0 noi isrslaxs ovht eutnt ag | 
 Rotistuges éAf ot 'efadbtvibet bec to ysatwr ert eo 
to s6Tun en? TmeQUS.<27 Nnkeeeosn af 3! ered 6 ‘ 
redmun rit ip Sikmitze Ae NG ise rese vil sudos erent 
Gevsm (All s78w aagtiy «nee fom Jud Inezerg ete tes 

neee don Sud duh s) ert 2S: ye@tis useq ort ot ineeengy alamins 
Squtt? shi ni emii- sane 2e& nese os 75 vert  teSwort edt 

to ofde7 er ya ct gnty gtd iamtve baasansnt 2f° yt Peigup eiat |) 
ee? 79296 VedL dette ot AWonm ef Gains cepsm to sedmuer ene , 
jaf? geet 46 cactygn <ahtrad) | nt - \o. AG) Berl aQree an | 
Tol ary ahs #4 (N'A beawuigete: vi (oeugaedam on) 

107 ome ett Nitgloa 4% berutobos) vl Maumsadue el subweb . 
ot), eslinnt tes. i AN nine an Oe ie) [Abs moett. ,eeéesto, ted aa 
ter Sod Mosmi? te driepeng 2! eubtwiba! bedum +e: tela! i: 

ats entelGo ene , (Ai\m2 of WISI eit prtbbs ve nent nenutaa |p 
(AAW eh Hol tmms eiaT (AIM Yo otamt tes boodt (eat iumrrKem 
- Pha saith mefomty HOMeniT tes storw. oii ot yao ater 
eet isaps nit teeth 789 eM nfems9 eA 30 tem mt) Be 


¥ 
a 


i 
oe 


7 a 


- 
_ 


ee 

_ 
7 
: i 
: 


- . | 
— | . | ~ 
or oan Oe ~ > 4 - re 


o2 


The equation for estimating the current size of the 
population takes the form of the simple Peterson estimate 


employed by Lincoln in 1930 (Seber, 1973): 


(14) we = Mx = ) 


lee ei Seeds 


Here, the number of individuals in the i’ th sample, Sik), is 
inflated by the proportion of marked adults in the 
population which are captured in that sample. 

The probability of capture is simply estimated by 
proportion of the population which were unable to evade 


capture at that time: 


(15) Pu = {ou = oat ) 
Ne TAte 
ee Ke= = (2-0) © 


The equation determining the number of young, which survive 
the ‘r’ years, to emerge as mature breeding adults at time k 
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The denominator of (16) counts all those individuals which 
were captured in and prior to the kK’ th sample and were 
subsequently recaptured. The summation in the numerator 
tallies the individuals which were in the (k-r}’ th sample as 
young (and are thus ‘r’ year old adults at time K) and were 
next seen in some K or later sample. Thus the expression 
represents the proportion of individuals subsequent ly 
recaptured that are ‘r’ years old at time kK. This proportion 
of the marked animals, Mik), provides the desired estimate. 
With the result of (16) the probability of the young 
surviving the ‘r’ years to adulthood is easily estimated by 
the proportion of the young marked at time kK which are in 


the population at time (ktr): 
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(kK+1) is the ratio of appropriate marked individuals: 
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This is simply the proportion of marked adults at time (k+1) 
which have survived from time k. However the young that are 


just becoming 'r’ year old adults at time (k+1) must be 
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discounted as this cohort was not part of the.marked adult 


Cessna fe time Ke 


The estimates for the number of unmarked adults and the 
young in the population both employ the Peterson type 


estimation equation: 
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Finally, repeating equation (12), the number of 
recruits (unmarked) to the population from time Kk to (k+1) 


is estimated by: 
(42) Bx = Miers a Ei, (Nic See) ) 


Forte ary c(t et) 


N(k+1,0) is adjusted for the number of unmarked adults 
already in the population at time (K+1) which have survived 


to time k. The remaining quantity is the estimated number of 
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Map Showing the Location of the 


Figure l. 
Anderson River Delta 
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Figure 2. sap of the Anderson River Delta 


Il. Data Analysis 


A. Background 

The brant is a migratory goose, breeding in circumpolar 
regions and wintering further south along the sea coast. In 
summer they inhabit loose colonies near salt water wherever 
suitable habitat is found. I recognize three subspecies of 
brant in the world, European brant (8B. b. bernicla), 
Pmldteive brant (5..0, mrota), and Paciviic prant be.0. 
nigricans), the latter two occurring in North America 
fearry, 19607). 

Atlantic brant, found along the eastern seaboard in 
winter from Massachusetts to South Carolina, move north in 
the spring to breed in the central and eastern arctic. 
Pacific brant, on the other hand, winter along the west 
E0asu Trom British Columoi1a to the Baja California 
Peninsula. They breed along the arctic coast from western 
plaska to the central arctic and. the high arctic 1s lands: 
Because pairing takes place during the southwards migration 
and on into the winter, and copulation occurs during the 
Hontnwards midration, there ais little interpneeding between 
Atlantic and Pacific brant where their summer ranges overlap 
in the central arctic (Barry, 1967; Bellrose, 1976; Palmer, 
1S Ie 

The data to be analyzed were obtained from a 
subpopulation (colony) situated in the Anderson River Delta 


on the Beaufort Sea coast in the Canadian western arctic 
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(Fig. 1). This banding program was conducted by Dr. 7. W. 
Barry of the Canadian Wildlife Service, with the assistance 
Sto coe cauthor., 

The Anderson River, emptying into the Beaufort Sea 150 
Kilometers east of the great Mackenzie River Delta, has 
formed a rich tidal delta covering about 65 square 
Kilometers of shallow channe1s and low flat islands 
(Fig. 2). The brant arrive there in late May. After the snow 
clears and the ground dries, they begin to establish 
territories and nest sites on the lower islands of the outer 
delta. Non breeding brant, yearlings and two-year-olds, do 
not ordinarily use the delta but spend the summer elsewhere 
in the arctic. After the-hatch in early July family groups 
gradually begin to congregate into a number of large flocks 
scattered throughout the delta and by the end of the month 
the adults have begun to molt. These flocks usuaily occupy 
the traditional molting sites of Kettle Hole Point, Brant 
felanc, Study Area island, Fox Den Island, wvwuncus sland, 
ieqancles 5 land. Loon Channel. Midd lecCnec sand sonearo in 
Creek along the brackish tide zone of the outer delta. 
Within three or four weeks the adults are flying again, 
somewhat synchronous with the first flight of the young in 
imescagly part of August. In Jate AuguS_wanclean ly) september 
the brant begin moving westwards along the coast, mingling 
with other groups as they work their way towards their 
common staging area in Izembek Bay, Alaska. From there al! 


migrate south to their coastal haunts for the winter 
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(Bellrose, 1976; Palmer, 1976). 


B. Methods 

The timing of the annual molt to coincide with the 
initial flightless period of the young provides an excellent 
opportunity to capture both the young and adults for 
banding. The brant’s propensity for flocking during the molt 
somewhat simplifies the logistics of capture but does not 
make it altogether easy. 

Banding drives were undertaken during.the last week of 
July and the first week in August when the young geese were 
large enough to retain a metal leg band and strong enough to 
withstand the stress of the banding. Brant were herded into 
circular banding pens, 10 to 15 meters in diameter, 
constructed of fine meshed fish net (herring seine) hung 
from poles 1.5 meters high and staked down to the ground. 
Two 50 meter long leads funneled the geese into the pen. By 
hazing about twenty birds at a time into a small! catch pen, 
within the main pen, banding personnel could easily capture 
and mark the birds with a numbered band placed around the 
lower right leg. The age class, either young or adult, and 
sex of each bird banded was recorded with its band number. 

When a previously banded bird was encountered in these 
Operations its band number, age class (adult), and sex were 
recorded. Often the recaptured band was so worn as to be 
partially or completely illegible; these were removed and 


replaced, and the new number cross-referenced with the old 
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one in the records. Most of these illegible bands could be 
read after being chemically etched, a process that removed 
the top layers of the metal, exposing the indentations of 
the numbers beneath the surface. 

During the twenty-one years from 1959 to 1979, as many 
as possible of the nine traditional molting flocks were 
sampled. There were, however, variations in the methods of 
herding. Between 1959 and 1973 a small motor boat and 
several people were used to herd and manoeuvre the brant 
into the pens. Wind, tide, the many shallow channels and 
hidden mudbars, combined with the rather stubborn nature of 
the brant, complicated the process of herding by boat. These 
factors limited the choice of banding areas sampled to 
Kettle Hole Point, Middle Creek, Shearpin Creek and Juncus 
Island. Since 1973 helicopter aircraft have been used to 
locate the geese, transport people, and to help herd. Within 
the limitations of funds and time this allowed a greater 
number of the brant to be sampled, especially from the 
me ling. fiocks.that were previously imaccessibje-=by boat. 

In 1968 and again in 1971 no banding operations were 
attempted. Furthermore, banding was not possible in 1974 
when all nesting attempts of brant failed due to predation 
by foxes (Alopex lagopus and Vulpes vulpes). As a result of 
thewtoxes’ activitiesSp»ciume acgul b.brantylentwinerdel tamonea 
imme emoLigmignd tion we hOrsthese, three. missing yearns. | have 
generated data conforming to the general trend of other 


years. For the purpose of testing my estimation equations, 
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the real banding data and the generated data will be 
sufficient. The problem of missing data in capture-recapture 
experiments is very real, and is open to further 
exploration. 

Over the years several different types of bands were 
used in an attempt to reduce the loss of bands from salt 
water corrosion. Until about 1967 soft aluminum bands were 
in use, which proved to have a very low resistance to 
corrosion, often becoming completely illegible within two 
years. In many cases the band was found to be so weakened 
tnaowe ftewase readyetomral loft theslegqeaSincer(96¢75 butywith 
the exception of 1973 when aluminum was used again, the 
bands were of harder alloys (monel and nickalloy). These 
have been more successful. 

Each spring at the Anderson Delta from 1959 to 1965 and 
from 1972 to 1979 a portion of the brant’s nesting habitat 
was searched for nests, the clutch size of each one being 
recorded. These data are used to estimate the maximum number 
of young produced by the Anderson River brant and provide a 


check on the results of the analysis. 


C. Results 

Inthe estimation equat vonsi tiromiPanmtil tinal anes used 
in this section r is given the value of three, the number of 
years required for the young brant to mature. marihestels, 11. 
and 111 summarize the raw capture-recapture data collected 


from 1959 to 1979. The distribution of the recaptured adults 
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which were originally banded as adults is given in Table |] 
for eden worethewbanding years ‘according tosthesttime ‘of “last 
capture. Table II shows the distribution of recaptured 
adults originally banded as young. These two tables, in 
effect, show the number of sampled individuals in each of 
the capture-history groups. Certain row and column sums used 
in the estimation equations are given in these tables as 
well. 

The number of the banded three-year-old adults just 
returning to the population is shown in Table III. These 
birds may have been more than three years old when they were 
finally recaptured, but they were all new mature adults at 
the appropiate time. Individuals of this sort which are 
recaptured when they are three years old appear along the 
upper diagonal, while those birds which are recaptured later 
in life appear along a diagonal deeper within the table. 

The only other remaining information required for the 
parameter estimates is the number of unmarked adults, 
S(K,o), and the number of young of the year, S(k,y), 
captured in each sample. These data are presented in Table 
IV, along with a condensation of the information contained 
ireiab lese a oli ands biteoresentecdians the tonm or ithe sums 
which appear directly in the estimation equations. Applying 
eqvicge Omseialc) throughs W000 domaine) ana lives SanmpParte i: 


with r set to three, results in the parameter estimates 


given in Table V. 
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In summarizing the more important and useful results, | 
have deleted the first estimate for each parameter because 
Enese were? Gale] ated With Yinsutt icient tdatam thatiis psthere 
were very few capture-recapture occurrences in the first 
several years of the study. The estimated brant population, 
N(K), at Anderson River had a mean of 4349 birds per year 
(ineluding young of the year, and both the marked and 
unmarked adults), ranging from 2801 to 8134 individuals. The 
marked adult portion of this population, M(k), generally 
numbered 814, while the mean number of unmarked adults and 
young was 1578 and 1957 respectively per year. 

The adult brant appear to have a high probability of 
‘surviving from year to year (mean A(k) = 0.70). Their 
Ghancesmot rsuUrvViving efor threetyeane) 1p = el) s4)mis 
significantly high as compared to the probability of the 
young surviving for three years (mean Y(kK) = .15). Unmarked 
brant were recruited into the population at a mean rate of 
808 birds per year. 

The clutch-size data collected at Anderson River is 
Slimmaritzed: imel able Vil J The ports onMot Mthne nesting@habitat 
surveyed each year is used to estimate the maximum number of 
young that could be produced, generally about 5600. However , 
the 5600 young is a gross over-estimate of the number of 
PaGse) ac tual yibavan lab le. for sbancangyy since tpreda tors: yiare 
many of the eggs and young long before banding commences. 
Moreover, the estimate of the percent of habitat covered is 


rough and serves only as a check on the analysis results. 
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Table VII shows that from 3% to 33% of the bands of 
recaptured birds had to be replaced each year, with a total 
of 230 replacements from 1962 to 1979. Eleven bands (5%) 
were so corroded by salt water that chemical etching failed 
to bring out the numbers. The life of the softer aluminum 
bands which had to be replaced is between three and seven 
years. Only (3%) of the harder bands used in the later years 
had to be replaced. 

Banding locations varied considerably, the most 
consistently sampled sites being Kettle Hole Point and 
Middle Creek. The remaining 7 saw more intermittent use. The 
distribution of the banding drives made each year is given 
me table Vill, 

The individual brant at Anderson River have a 
propensity to nest and molt within the same areas of the 
delta year after year. To illustrate this phenomenon I have 
tabled all of the recapture occurrences by location, 
according to the location of previous capture (Table IA). 
Many of the nine banding locations are in close proximity to 
one another so I have recognized four nearly autonomous 
brood rearing areas within the delta. Examination of this 
data suggests that intermixing between sites is low, 57% of 
the capture-recapture occurrences take place in the same 
avease) i reating (able 17 asa four by (eur scontingeney table 
and applying a Chi-Square test for independence indicates 
fhat recapture location is very dependent on the site of 


previousmcapture (Chi-square = 9/3/30, dat sa 
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The source for the following tables is the unpublished 


field notebooks of Dr. T. W. Barry and the author. 
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Capture 
Year j 
1959 
1960 
1961 
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1963 
1964 
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Jabie liv Bandingevistribution of S(1, jy} 


Recapture Year ji 
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Year j 1960 196 1 1962 1963 1964 1965 1966 
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Table II: 


Capture 
Year j 
1959 
1960 
196 1 
1962 
1963 
1964 
1965 
1966 
1967 
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1969 
1970 
1971 
1972 
1973 
1974 
1975 
1976 
1977 
1978 


Cont inued 
1974 1975 
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0 0 
0 0 
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0 0 
1 0 
0 0 
2 1 
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habte 111: Continued 


Year Originally Banded as Young, (i-3) 
Re- 
Capture 


Year j 1965 


1968 5 1966 

1969 7 6 1967 

1970 2 g 1968 

1971 0 4 g 10 1969 

1972 3 14 g 3 1970 
1973 0 2 11 3 3 2 
1974 0 0 0 5 O 3 
1975 0 ! 5 0 6 0 
1976 0 0 2 i 0 0 
1977 0 0 0 
1978 0 0 3 0 2 
1979 0 0 2 2 0 0 
Sum ate) (15) (54) Cen (14) (7) 
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Table III: Continued 


Year Originally Banded as Young, (i-3) 
Re- 
Capture 


Year i 1971 


1974 13 1972 

1975 11 14 1973 

1976 3 6 5 1974 

1977 2 5 1 9 1975 

1978 4 5 7 2 8 1876 
19783 2 3 2 9 - = 


Sum (35) (33) (15) (30) eer’ (9) 
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Table IV: Summary of Banding Results 


Banding Year 
Samp le 
Observ. 1959 1960 1961 1962 1963 1964 1965 


Sikyo) 2 105 19 264 471 158 825 
Sikey) 22 246 DZ 334 447 170 491 
Sm(K ) 0 0 Q ho 54 26 74 
S(k) 29 sou 78 613 972 354 1090 
* 1 3 26 33 i418 123 
x 7 ) 5 18 28 
EK 0 ‘i 0 5 
KKK 1 30 10 61 106 84 136 
Z(k) 24 19 18 95 95 
R(k) 1 30 10 56 $3 69 119 
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Table IV: Continued 


Banding Year 
Samp le 
Observ. 1966 1967 1968 1969 1970 1971 1972 


S(k,o) 401 319 357. «185 40 157-269 
S(k,y) 236 555 254 236 56 289 458 
Sm(k) 89 135 134-160 115 ie = tee 
S(k) 726 1009 745 #581 211 589 = gg 
x 170-265 369-406 429 326 = -:263 
** 32 69 85 45 ot 40 68 
EK 13 1 17 15 54 a4 14 
xxex 230 238 RIT: 40 104 242 
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Table IV: Cont inued 


Banding Year 


Samp le 
Observ. 1973 1974 1975 1976 1977 1978 1979 


Sikeo) =.256 251 403 119 67 406 400 
SikK,V) 367 396 Hehe} 308 250 928 awit 
Sm(k ) ie 208 291 22) 115 286 33 
S(k) 795 555 1399 648 432 1620 1408 
* 333 400 Dar Gish 364 141 0 
** 48 45 62 4 
HK 7 35 33 is cia 32 9 
RK DAES 174 289 128 63 190 
Z(k) 285 sists 225) 310 
R(k) 260 144 pda ope 
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Table V: Parameter Estimates 


Banding Year 


Para- 

meters 1959 1960 1961 1962 
M(K) 69 110 
N(k) BiG 4482 
P(k) O.10 0.14 
: 0 
¥(k) 0.00 0.06 0.00 0.05 
A(k) 0.38 
N(k,o) 188 1930 
B(k) i Ou ae 
N(k,y) 516 2442 
= Dore 2040 
ae 0.66 0.54 
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Table V: Continued 


Banding Year 


meters 1966 1967 1968 1969 1970 1971 1972 


M(k) 404 619 974 1164-2031 1378 574 
N(k) 3292 4624 5415 4228 3727 5675 3072 
P(k) Tee Bee Dra 0.42) SoeOGee Oden 1 0no9 
* 23 28 40 64 eG OC pee? 

¥(k) ee Wee OOCIo Ra OREOT® 0.18 Die? eee: 
A(k) DAio 7) Oho Ones 4 0.62 0.36 0.75 
N(k,O0) 1818 1462 2595 1346 707 1513 924 
B(k) 422 1456 0 673 1101 438 464 
N(k,y) 1070 2543 1846 1718 989 2784 1574 
x Spe pais 3569 2510 2738 2891 1498 
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Table V: Continued 


Banding Year 


meters 1973 1974 1975 176 URS ef 1978 1979 


Mik) 641 1340 888 1107 
N(k) 2964 5506 4268 3245 
P(k) C27 aewte Oey 0.20 
* 10 83 57 31 
Y(k) 0.85 

A(k) Os) a OEIC 

N(k,o) 954 1617 1230 596 
B(k) 639 517 0 

N(k,y) 1362 2549 2150 1542 
> pei 2957 2 ahs iW oc 
4% 0.46 10.46 0.50 0.48 
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Table VI: 


Samp le 
Year 
1959 
1960 
1961 
1962 
1963 
1964 
1965 
1972 
1973 
1974 
1975 
1976 
1977 
1978 
1979 
1980 


Clutch esizes. 


Number 
of Nests 


221 


Anderson River Delta 


Number 
of Eggs 
692 
1518 


Habitat 
Area 
Sampled (%) 
10 

as 


Oe 


Estimated 
Number 
of Young 
6920 
6072 
6970 
6330 
378 
4170 
6003 
3740 
5820 
3480 
387 
5100 
8360 
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Table VII: Distribution of Replaced Bands 


Samp le Rep laced Unetchab le No. of Recaptures 
Year Bands Bands Recaptures Replaced % 
1959 

1960 0 0 0 

1961 0 0 0 

1962 5 0 HS ce 

1963 13 2 54 24 

1964 3 { 26 12 

1965 14 2 74 19 

1966 ee 1 89 0) 

1967 26 0 125 19 

1969 18 1 160 11 

1970 3 2 ake 3 

1972 28 0 ro 14 

1973 ¢ 0 fig 5 

1975 17 0 Pe) | 6 

1976 13 0 224 6 

1977 13 0 is ia 

1978 36 2 286 13 

1979 24 0 eter is 
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Taolies VIL Dirstrabution of Banding Locations, 
Anderson River Delta 
Group 1 Group 2 
Kettle 
Samp je Hole Brant Study Area Fox Den 
Year Point Island Island Island 
1959 > 
1960 * = 
1961 * 
1962 * * 
1963 - 
1964 
1965 
1966 = es 
1967 He 
1969 * 
1970 = 
1972 * 
1973 a8 
19975 *: ’ 
1976 - a 
1977 * 
1978 * is 
1979 * i 
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Maple Viil> Continued 


Group 3 Group 4 
Samp le Juncus Triangle Loon Middie Shearpin 
Year Island Island Channel Creek Creek 
1959 
1960 
196 1 
1962 
1963 = 
1964 = * 
1965 - 
1966 * 
1967 < = * 
1969 = 
1970 
1972 * 
1973 = 
Weites 
1976 * : 
1977 = 
1978 1 “ 


19793 * i" 


Table IX: 


Capture 


Area 
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Creek 


Use of Molting Sites by Brant, 


Anderson River Delta 


Recapture Area 


Kettle 

Hole FoxDen Triangle 
Point Island Island 
365 202 48 

87 162 62 
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92 68 24 
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Dv Discuss 10n 

This analysis of the banding data of the Anderson River 
brant provides reasonable parameter estimates. My results 
indicate that the number of adult brant (N(kK,o) + M(k)) 
breeding in the delta ranged from a low of 1513 to a high of 
4470, with a mean value of 2390. This interval encompasses 
the estimated population size of 3000 for this colony 
(Barry, 1967). Barry's figure was derived from direct 
observation, clutch-sizing operations, and aerial surveys 
during a period of years overlapping with this study. 

My mean estimated yearly production of young (Nik,y)) 
is 1957, substantially less than the possible number of 
young produced from eggs at Anderson River each year (Table 
V). However, this colony has a high egg and chick loss due 
to predation and natural mortality prior to banding season 
(pre-fledged stage) (Martin and Barry, 1973; Barry and Barry, 
in press). A mean production of 5600 eggs is reduced by 23% 
to 27% (these percentages are the result of clutch size 
studies on brant in other areas), (Bellrose, 1976; Palmer, 
1976: Pacific Waterfowl Flyway Council, 1978) to about 4000 
Suiccesstully ‘hatched and dried chicks. “Furthermore, ‘during 
the next three weeks prior to fledging, these survivors are 
further reduced in number by about 30% to 35% (Barry, pers. 
comm.: Bellrose, 1976), leaving only about 2600 young 
available for banding. Further assessment of my estimated 
number of young is made by examining the year ly percentage 


of young produced by the Anderson River brant, derived from 
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the estimates for N(K,y) and N(k), a quantity ranging from 
26% to 55% (mean of 45%). This is substantially above the 
10% to 40% range found for the total Pacific population 
(Pacific Waterfow! Flyway Council, 1978). However, local 
conditions at Anderson River, such as early snow melt 
Relative ito ithebrant ss ‘anrival.,. longer daylight, aless 
precipitation, and generally more clement weather, could 
favour a higher percentage of young to be produced there 
than at other colonies, e.g. Yukon-Kuskoquim delta (Barry, 
pers. comm.; 1967). Furthermore, estimates of adult young 
hat tos are «made .-in the winter, after the fall) migration and 
after natural mortality and hunting harvests have taken 
their toll (Pacific Waterfowl Flyway Council, 1978). 

Analyzing band returns of brant breeding in the Yukon 
Delta, Alaska, Hansen and Nelson (1957) arrived at a 45.4% 
mortality rate for young brant in their first year. They 
also reported an average mortality rate of 32.2% for each 
succeding year. In terms of survival, these figures are 
equa Valent sto. 25914 of the young sunviving athe mtnnee years 
to adulthood. Their result is substantially higher than the 
15% mean survival rate resulting from my analysis. The 
mortality figures of Hansen and Nelson are, however, derived 
from a simple deterministic life table method developed by 
Bellrose and Chase (1950) who express the mortality rate as 
the proportion of banded birds which died during the year as 
determined by bands returned by hunters and bands found on 


dead birds. The obvious deficiencies of this method are that 


. i 


blues , 1eMi. sew « aL ened bre eo 34) tat ena Pi 
ayedt besubo1g ad oh ee 72. 90s insoteq rerieid 6 ‘quovs? “ 
wangse' sitéeb mrueowsun= Roney ele ,zatnefod write sae = 
gaulay T1Gbs. To cesantsae SromAantasF . tact AOS. 279Q 
baie netisrtgtm ifet ara nats oéiotiw set of ebsm ets 2otist 
NSAes aver asesv 1er patter ois vit 'so7om Lesten n91ts ie 
EST9t) A uBd veel 2 lwot regal ot Tee Sd het ited) 
nov’ ent ot pat beard 746d do anise 1. brea pareytend ve 
Shi Sh B icles oavints (Teel) apelat pre osensh., shells «siied - AP 
WSnT.-. qaouvr T2q7tt. speci nh tqsadegmvoy jet ets “iristsom 
Aube noH-BO.52 to ets7 viii strom Soe evs ne -belneges oats 
ane 29qupth aeons wievivque To amie? ol , aay gritbesoue a 
askey seins ets pervivavue Onboy artz 76 Sf 12S of inal evtupe as 
adi os? sendin yl lst tnsisdbe ef tfueer + slid Leslie cil 


bevieaet ..1evewon 945) noel si bas nsensH +o wen) a te 
ye hegolisveb boniam aides ime? ofteintareieab signe 6 
a6 9f87 Yi! si som ont) 2eetaxs ortw. (969) ) eget bra. 
as: They, Bd?) Ooi 7ub bsie ndtdie : 


- 
7 


og 


not all band recoveries are reported by hunters and 
certainly not all dead birds with bands are found: that is, 
dead banded birds are found with a lower probability than 
the dead unDanded birds, since there are relativly few 
MarKecebirds im thewoopulation®’ This®bias restsitseinean 
Uncgerestimatvesore the acttal mortality rate. eies. wan 
everestimate orythe@sirvivalWrate: 

For the adult survival probability a more reliable 
estimate than that provided by Hansen and Nelson is given by 
Boyd (in Le Cren and Holdgate, 1962) and lends credence to 
my estimated survival rate of 0.70. Boyd applied an 
estimation technique, developed by Haldane (described in 
Seber, 1973), to the banding data of Hansen and Nelson and 
found that dasie fy aver Dbrant had@a constant survivel! pate 
of 0.85. Using banding data collected in Spitzbergen, Boyd 
also calculated that the closely related Atlantic brant 
SUrVYiVeGewrin a probabi Tityor 0983" “halddne s@siochastiic 
esuimattom precedure Ss morewcred ib eV inant tna tet Se ltrose 
and Chase since he at least included in his analysis a 
Gomstant probability that a band is found and reported. 
However, Haldane’s method is not completely realistic, since 
tats probabi l7ty coulic *chamge=w7 th "hunting pressure anc 
banding effort; for the purpose of comparison though, Boyd's 
estimaces wi) be suTT TCreltt: 

ihe remaining Pepulation=pdhrdameverc arem Naki K= aay yak 
the number of three year old adults, and Bik), the number of 


recruits. There are no estimates of these parameters from 
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other sources with which to draw comparisons or make 
assessments. 

several aspects of the assumptions must be examined in 
order to ascertain the validity of applying the estimation 
equations to sthiseparticulareset.of.data. —Asumotedsin the 
results, the Anderson River brant are likely to use the same 
specific molting areas (these being Kettle Hole Point, Fox 
Den Island, Triangle Island and Middle Creek) within the 
delta from season to season. It is evident then that the 
marked geese, released after each season’s banding drives, 
do not thoroughly mix with the remainder of the population 
by the following year or, for that matter, possibly ever. In 
addition, there were different sampling techniques employed 
during the course of the study (boat or helicopter). The 
concentration of sampling in only a few areas versus the 
later general sampling by helicopter causes problems. It is 
very probable that not all of the Anderson River brant have 
the same probability of being captured in any one season and 
that the sampling is, in effect, non-random. Brant that are 
banded on Fox Den Island or Triangle Island, for example, 
might have little or no chance of being recaptured if that 
site is not resampled in subsequent years, even though the 
birds may be present in the delta. 

The effect of this non-random sampling on the estimates 
is twofold. Intuitively, the estimates of the survival rates 
of both young and adults will be biased on the iow side if 


enough individuals are involved in one of the infrequent ly 
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used banding sites. Since the survival rates are basically 
estimated by the proportion of marked birds that survive to 
and are captured in the next sampling period (third next for 
the young), it is clear that this estimation procedure wil] 
interpret the event ‘a greater number of marked birds not 
being recaptured due to non-random sampling’ as the event 
that ‘more of these marked birds are dying before they are 
recaptured’. As a result of this misinterpretation a lower 
estimate of the probability of surviving is produced. 
Secondly, the estimates of the population size will] be 
inflated since more marked geese could potentially be 
captured if the sampling were random; the estimation 
procedure translating (see equation (14)) the event ‘reduced 
proportion of marked in the sample’ as the event that ‘more 
geese are in the population’. This results in an estimate 
that is biased on the high side. 

Another source of error in applying these estimation 
results to this specific population is that assumption 5 is 
violated: geese do not lose their bands during the course of 
the study. The rough estimate for the rate of band 
replacement and band life represents minimum values for 
these statistics and suggests that the band loss is very 
high from year to year. When the softer aluminum bands were 
used, 11% to 33% of the recaptured aluminum bands had to be 
replaced each year. That this problem occurs and complicates 
the data analysis is very certain, the ramifications of such 


high loss being inflated population estimates and deflated 
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survival rates caused by a reduced proportion of marked 
individuals in the population. The reasoning of this 
argument follows exactly that given in the previous 
discussion, the effect of band loss and non-random sampling 
are identical. 

A further problem affecting the analysis of this brant 
population is that on the spring migration the recently 
matured males follow their females to the colony where she 
originally hatched, seldom returning to the colony where the 
male was hatched. Thus at the Anderson River there is a loss 
of banded males and a corresponding influx of unbanded males 
following the returning local Females. ihiseresulits ind 
dilution of the supply of marked birds and hence biased 
estimates. examination of band returns has shown that such 
an exchange exists between the brant breeding in Alaska and 
at Anderson River (Barry. 1959 to 1979) and is also Known to 
occur with populations of snow. geese (Chen hyperborea) 
(Cooke et al., 1975; Finney and Cooke, 1978). 

I will also mention in passing (and without any direct 
evidence) that there is a distinct possibility that the 
tacit assumption of brant always breeding at age three may 
not be valid. While preparing the banding data for Tables I, 
Peeanceel lle ienoticed that dumingmtice les teurnree los nOuG 
years of the banding program a significantly higher number 
of one- and two-year-old geese were being recaptured for the 
first time, but that overall, there was a tendency for the 


young. brant to be first recaptured five. six, or more years 


: 
inetd gin to 2taylsnk oat Mitioat ts mest a eanhaest a 
vitqsesy ant err eer adi no tert ai patistug 

ette Sagrtuewteips, 90) ot pelea ated wot ig? aetem. be ut sm 
sit gtarwoviielioe) S44 ot Qainotg mdb se, bedtgert a 
2am! 6 ef Gast. hevrh meevetria #at te eunT Jberisted saw 3 = 7 
agtam behnecitu Vo subi) peiynGqa co sp bas esiom babnad, 3 ° i 
5 mt siauees eta! aatanet pool guiorwites 202 articles 
bezeid sans brah ree ra, tates > ylaque ans > net tulip. 

Aowe Tartt fworte ean arn &> Pate tc nohlantaacs eetemt te iv 
one edenre at otpsbad irsne of) /neowied 2/ai re egnenoxet 
ot Awonm dais ai Bins QETET sad Reel. vars 1evtH mozsgheA 1B) 
| Rye oleh qand) SessoJwone Io. snot isiwan mi tw 4 990. 

(@%er sano? come Vefart? 7¢CTet Aa Ie a7009).° _ 

logestS yng SuGlitw ome prised ei fortes oale Tiel 7 
ent dedi y*l Cidtgeo Tenfigib © et san: Der — 

vant Ser) SpA, 16 oni Gesrd: 2vew!ls tris, to cat hegre, 2 
J) ge?ae? 167 stab, gal ond enti ecinegend. 61 haw ‘ahtae 
wa? Of Seti? tes/ ant periye Jett beat tom i. 

as Cg) weer ¥! nysat? tate ‘% mA7go7, quad i To 
salad Mike seine! ied Lene 


later. 


63 


64 


Conclusions 


To briefly summarize, I have addressed the problem of 
estimating parameters of a non-homogeneous population 
consisting of two recognizable age groups, adult and young, 
where the latter require some specified constant number of 
years (r) in which to attain adult status. The estimation 
procedure derived for this type of population allows for one 
probability that the young survive this maturing phase and a 
cifferent probability that the adults survive the one year 
intervals thereafter. The estimates of population parameters 
are found by maximum likelihood techniques. 

In theory, the area of applicability of this estimation 
scheme are those populations not covered by Pollock’ s 
general analysis discussed earlier. This includes those 
species of birds, mammals, fish, et cetera, in which the 
subadult members may or may not be distinguishable from the 
adults but they at least spend the breeding season in areas 
away from the adults (that is, if marking takes place during 
the breeding season}. In practice, however, there are 
several factors which severely limit the extent to which 
thisedhd | ysisew)) |i oracould ibe used: 

Primarily, in thiseprocedure, more time required sro, 
the young to mature, that is, larger values of r. will 
require more samples to be taken. This is because the number 
of samples must be at least r before any calculations can be 


made at all. Furthermore, for any reasonable measure of 
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accuracy and reliability to be attained in the estimates, 
the marking program must be continued for a substantial 
period of time after the r’th sample, resulting in a 
requirement for a large body of capture-recapture data. 
Satisfying this requirement presents many problems with 
conflicting solutions when dealing with a species for which 
the young require a long maturation period; accuracy 
requirements must be balanced with the available study time 
and the urgency of the situation. 

Perhaps even more restricting in terms of the 
applicability of my analysis is that there are relatively 
few species which can be easily captured on a scale large 
enough to allow sufficiently sizeable numbers of individuals 
to be marked. There are even fewer species of this type 
whose populations follow the pattern required for this 
analysis. This estimation procedure is designed to be 
applied primarily to geese, about the only species meeting 
the specifications discussed above. However this designed 
application may be limited as well. For the forseeable 
future the era of long term observational investigations is 
drawing to a close as government funds for these types of 
studies are curtailed and even withdrawn completely in 
favour of investigations of shorter duration. 

ine results of applying this estimation procedure =1o 
the capture-recapture data of the Anderson River brant are 
sufficient to determine that the estimation equations are 


generating reasonable estimates that are near values 
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predicted through experience and from other sources. For use 
in management programs, however, the reliability of my 
results are not completely satisfying. This is, to a large 
extent, due to the lack of consistency in data collection 
procedures and to the presence of three years of missing 
data. Recommendations for further studies then would include 
a designed sampling method that would provide consistent 
sample sizes and would be adaptable to idiosyncrasies of the 
population in question. As well, effort should be expended 
to insure that tag loss is kept to a minimum. 

There are several inadequacies with the estimation 
procedure as | have derived it in Part I. No accounting has 
been given for the possibility that the survival 
probabilities are dependent on an individual’s capture 
history. For example the young brant that are captured and 
marked are released under a great deal of stress having 
Spentsup to’ four hours under cold iconditions: inmthe* pen, 
after being somewhat trampled by the flock as they are are 
hazed into the catch pen. On release the young often run 
away in many different directions, becoming widely scattered 
and separated from the parents and thus are easy prey for 
the glaucous gull (Larus hyperboreus) (Barry, pers. comm. ; 
Barry and Barry, in press). Obviously these young must have 
a significantly lower probability of surviving than the 
young which did not undergo this trial. 

Related to this factor or “a *caplure-nistory dependent 


survival rate is the very real possibility that capture 
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probability is not only time dependent but capture history 
dependent as well. A clear example of this is the 
white-fronted goose (Anser albifrons). When banding molting 
sub-adults in their traditional flocks, it is often very 
Cis teu tetol capturesmany it: thatesitelwasesamojedi im the 
previous year. On several occasions, a helicopter and a boat 
as well as several people on foot were all used to herd with 
the result that only a few tens of the 500 or so geese were 
actually forced into the pen, the rest disappearing, being 
led by the more experienced previously marked birds. 

Incorporation of these two options into the analysis 
is, in essence, measuring the effect that the observer (the 
bander) has on his experiment (the geese), an endeavour 
which can not succeed. Since the act of observing (taking 
samples) itself drastically affects what the experimenter 
will observe, it is uncertain as to whether the parameter 
estimates are representative of an artificial population 
that is created by the observer’s presence or whether the 
observations are made from the real population. In other 
words, what type of reality is being measured. This is the 
major problem confronting this sort of analysis, and is one 
that appears to be ignored by other researchers employing 
the Wek ener type of approach. 

Also neglected is any method for calculating variances 
and covariances of the estimates. Jolly (1965) produces 
variance estimates for parameters of his homogeneous 


population analysis by US#nG -4) truncated hhaviomiseries 


inal attr .e 
edt ni He hanse caw. orf 


rst a aa = ie oe 
Jeed 3 WAR a9iqoa! lst # . enaTenogo: Iensvee. 1 wa 
piste bearded wh eee “the 2 aaunen a ‘at goed heGoues ast 
grew eeaany Je co W9e ont Toe payed. wet ey noche vias 
eolsa (oni wwsaguaath fea adf. fq, seit ofa heat? t 
abatad ofdem Vi gueiverq beonsiisare@ e40n poe: se 

atava gra oni3 aint. 2Aetsae ous Saaqrt fo hat 1S4nqioanl . 
a02) Sev7Gado Su?) Jerid ioe ts eq! anhruaean Soeeas nm 
quevgstne ns.,'ea090 artti-taaninsexes ett no ean a. 


gniiangs) oni Vaeada. 16 Se il sorte hemoug ju#t ABD aan 


ive 


4sinamtrodes oft rirlvy ‘otaetta +i [nastent Steed lastone 


jdiemetea efit saith ben ee 2h ibd iow 2 Ph ,ee"7sedo ih 

noael.iqes Us! 3! bra ab ko Gy! leerete7ge9 Sip eetemt : 

ari ~“Griscw io soreeery 2. \ a acein adi yo Sefapya ef ff 3 

tito, foi dehiySet hgh aA, mt shen 646 anotsevreede 

ods. @) eff! peyuekem cAraaal ys riser Wa. one Sart e 27aW 

one 2’ hee efavisns to. 7 4o2 aefit anit Tus Tews: ‘tice | 

Qrigalame evotsieeest Jeri. qa Gerenyt ed of aasegqe 

‘Aewerays 10 Squ) wey? 

agarintivny pyedone Ml nai 0%, Sida “engi Hage 

wats cdl oie aa hd 


68 


approximation. Stokes (1980) and Pollock (in press) both use 
this same method of approximation in their analyses, and at 
the present time it seems to be the only means available for 
computing variances in this stochastic estimation setting 
The main criticism and objection to using and relying on 
these approximate variances is that they are asymptotic, 
requiring a very large sample size in order to return any 
accuracy. Furthermore, for small sample sizes the 
distribution of the maximum likelihood estimators are quite 
skewed (Robson, in Johnson and Smith, 1968) leading to 
inappropriate and misleading confidence bounds. 

Still open to investigation is the possibility that the 
estimation procedures may be improved by using data from 
other sources. One option is to combine the recovery type 
procedures of Brownie et al. (1978), which I discussed in 
the introduction, with the recapture method which 1 have 
developed in Part I. This would use the large body of 
recovery data collected for most species of North American 
geese and ducks. The estimation analysis could be modified 
to differentiate between natural mortality and that induced 
by hunting pressure, by introducing probabilities that young 
and adults survive the annual hunting barrage in the fall 
and winter. As well, the shot death of a banded goose 
definitely terminates development of the band’s associated 
capture history vector and so precludes any possibility that 
the band is ever recaptured. This certainty of information 


should have beneficial ramifications on the existing 
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Appendix: Derivation of Estimation Equations 


I will now continue with the derivation of the 


estimation equations following the derivation of equation 
(12) in Part 1. The final versions of these subsequent 
equations are given as (13) 


iy 


through (20) at the end of Part 


The number of unmarked adults at time kK has the 


Following likelihood function: 
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Setting AlIn L(N(k,o)) = 0 and using equation (12) for B(k) 
and B(K-1) and solving for N(k,o): 
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Thus the number of unmarked adults is estimated by the 
number of unmarked adults in the sample inflated by a factor 
fe AP Kas 

Using the estimate of the number of recruits, B(k), 
from equation (12) the estimate of A(k), the adult survival 


rate, given by equation (11) is revised as follows: 
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which reduces to the following after using equation (1): 
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This estimate finds the proportion of those marked adults at 
time K who have survived to time (k+1). The number of young 
that became adults at time (k+1) are removed from the 


numerator as they were not members of the marked population 
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From (8) the likelihood function of the young at time k 
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Like the unmarked adults, N(k,y) is the number of young in 
the kK’ th sample scaled up by an ammount 1/P(k). 
The likelihood of the number of young that are just 


becoming adults is: 
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Again, by setting the first difference of the logarithm to 


zero, one gets: 
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Several estimation equations will now be extracted from 


the likelihood function (8). By combining these with some 


previous equations an estimate for N(k+1,k) + N(kK+1,k,y) 
will be found. 
Firstly, estimates of those adults last captured at 


i ee! 


some previous time ‘ j’ and initially marked when adults, 
will be made. 

There are two cases to consider: 
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A description of the intuitive equivalence between the 
denominators is similar. The proportion of N(k-1,j) that 
survive to time k and and form the Ni(k,j,w) is A(kK-1). 
(1-P(kK)) is the probability that these survivors will not be 
captured, and so those that were captured (S(k,j,w)) must be 
removed from the left side. These uncaptured survivors die 
with probability (1-A(kK)) before time (kt+1). The survivors 
(N(K+1,j,v)) are therefore removed from the left side 
maintaining the equivalence. 

2) When j = (kK-1) there are two forms of the capture history 


w to consider : 
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Next,estimation equations for those adults initially 
marked as young are developed. As with the adults marked as 


adults there are two cases to consider. 
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a) For jee eee ay exceDt storethemeacemmt = miic-7i 
and w = t, the likelihood function is: 
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For large N(k,kK-1,y,w!, setting f\In L{N(k,k-1,y,w) to zero 


results in: 
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The intuitive description of the equations (A8), (AQ), (A10) 


and (Aii)} will follow that given for (A7) with few 
modifications. 

Since the right hand side of equations (A7) through 
(A11) are independent of the capture histories w, the 
numerators and denominators of each of these equations are 
summed separately over all capture histories w. 


This results in the following equations: 
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Since each = signt hand@s1tde tor A 12 threugh AiG ets 
independent of j all of the terms in the numerators, 
dependent on j, can be summed over the range of j for 
ath te ee Similar i yeral | mternomunie: ie 
denominators are summed over j, resulting in five equations. 
These five are reduced to one by adding the numerators 
together and adding the denominators together and equating 
Ene-result ‘to thesright ‘hand «sideYy Equation (615) 41s not 
Valid, however when j = (K-r) and w = t and so, this 
equation must be removed in the final analysis. Even though 
the denominators of (A13) and (A14) are the same, they must 
not be counted twice as they were derived from two different 
forms of capture histories. 

Manipulating (A12) to (A16) in thé manner described above 


will give: 


kez L 
Nec “1 Sg - Vi; + Sie pS Ren Day eee 


J=! j= { j- yo 


hat it y Se-iyy a Mary = Mi gene ryt + Sr- (Metis = Wk k-eyt 


yet 


Cen eee = a nanan ea ee ale 
Nj -Y sy - Tey cedtlne Ghee A mikace ony YE Sxy- 


Je j= ae je! 


Nusyy + Mrekevy ia Skey = Wee ke-ly Bice “rye + Skk- cyt +NMarnery t 


j= 


() Se) 
Cle Peek) Alen 


ES Nee gen (0 a ae 


anot fauiie avi? at pettus 
exgisnam oft nies « ot beouuny ate / 
got isups Ber sertaget sacishimbadh od) ‘priate bes 
166° 24 VerEI nerd nastes . Shia hiner sighs art or ity 
ale? be) beh oe ikenbh Soot a ogre -vevewor < 
Miguel? nave steylenk tant? ea? of bevomes sd teum not 
foun Yat acee eds asp (NAF bis, (etd r te anotentar 3 
irigag?? My Ow! mon? bavinah érew vets as. aohwi peznuos od 
.2gh7o3eid @sulgss, to amy 
ovetie bsdtisees vandbwert<ht. Chal of (Ona) qabvanughe 
ie Levig = 


~ pale ce god 3: weds allt.* vende + hs - lt vid oa) 


taser AY = aarty AZ yyetealh ~enelh - me + galt 1 
5 | 7 % a ; ; 
.< ans 4 te nah pane. = pata nl 2 - yas at” 
- . > 


~ gine : 


‘a yranalh+ tara ne “4a axel . ean = yravthe os Se 
| i | 


85 


The last three terms in the numerator and the denominator 


account for the case j = (k-r) and w = t from equation 
Ae One. 
Using the relation {1), the numerator reduces to 
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Sincer N(K-Wak-f Val) = ool hn Ker y. tl) J] 0 these anelst inl 
young and are assumed to be out of the population). 
The denominator is similarly rearranged and after 


adding and subtracting N(k+1,kK) + N(k+1,k,y} becomes: 
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And so, the equation relating the ratio of probabilities to 


the ratio of the corresponding expected frequencies becomes: 
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The numerator of (A17) can be described as follows. The 
first term is the number of marked just after the (k-1)’ th 
sample. The second term are those marked adults at time kK 
except for those ‘r’ year old adults who have just arrived, 
and therefore were not among the marked adults at time 
(K-1). The difference then is the number of marked adults 
that did not survive to time k. Divided by an appropriate 
quantity, this will estimate the probability of dying before 
time k, (1-A(k)). 

In the denominator, the first term are those adults who have 
survived from time (kK-1) corrected for the ‘r’ year old 
adults just arriving. From these are removed ineee marked 
individuals who were in the k’ th sample, again correcting 
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denominator now describes a certain group of marked adults 
who were not captured at time k. After accounting for the 
‘r’ year old adults at time k and the marked adults that 
were captured at time k, the third term removes those marked 
Ages taal SUrVIVea toutime (Ki), ihe Quamerty then lent 
in the denominator divided by an appropiate term estimates 
the compound probability in the right hand side; that is 
‘survives to time (k-1)’ and ‘is not captured at time k’ and 
‘does not survive to time (k+1)’. | 
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right hand side of (A17) becomes: 
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Using this form of the right hand side and rearranging, 


(A17) becomes: 
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This is then simplified further, by substituting with the 


estimates given in equations (10), (A2), (A5) and (A6), to: 


(A18) Mice = ae Nrsiry = Ax ( om Gu) 
eee ie eels 
This simple equation equates the estimated survival rate to 


the estimated proportion of the marked animals released from 


the k’ th sample that have survived to time (k+1). 
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The last step in finding maximum likelihood estimates 
will be to establish recursive relations for the Nlciy 3 Ng 
and the N(i,j,y,w). From these an estimate for M(i) will be 
extracted which will be applied in the necessary equations 
developed earlier on. 

Wit K@sev torel+i) tam 8 ie the kel thoce=functionrof 
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gives, after summing over the capture histories w: 
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Setting Kk = ) and rearranging, equation (A12) becomes: 
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Then using (A19) in the right hand side and multiplying both 


sides by 
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to both sides results in: 
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Multiplying both sides by 
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the numerator of the left side counts all those adults at 
time i last captured at time j who did not survive to time 
Mitim. analogous, to, the numerator of the right side, the 
probability of dying before time (i+1). The S(i,j) must be 
removed as they are no longer members of the N(i,j), after 
capture they become members of N(i+1,i). The denominator, 
meanwhile removes from the N(i,j) all those who were 
subsequent ly recaptured at some time up to time |. This 
leaves the number of the N(i,j)’s which are never 


recaptured. Again, divided by N(i,j), this gives an estimate 
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of X(i), the probability of not being captured after time i. 


Again, with k=(1+1) (8) contains the likelihood 
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For large N(1+1,j,y,v), equating the first difference of the 


logarithm of this function to zero gives: 
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These equations (A24), (A25) and (426) will be used 
successively to develop the recursive relation for the 
Niet Gy Ww). 

In equation (A10), let k=1 and sum over all capture 
histories w, except for w=t when j=(1-r) which must be 
excluded. When j=(1+1-r) and w=t the equation is present but 


N(1+1,1+i-r,y,t) is never given as the last term in the 
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denominator in the left side (since all of the terms in the 


numerator in this case are zero) and so it must be excluded 


as well. 


Thus (A10) may be written as: 
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Let K=1 in (A10) as above but when summing over capture 
histories exclude the case w=t when j=(1-1-r) as well. By 
following the same sequence of steps that lead to (A27) and 


using (A25) in place of (A24) we have: 
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Intuitively this equation for adults marked as young is the 
same as that derived earlier for adults marked as adults 
(A22) with one change. The final sum in the right hand side 
of (A31) removes all of those individuals who were still 
young at time i and as a result could not be in the i’th 
sampling. 
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The denominator of (A36) counts all of those individuals 
which, captured in or prior to, the k’ th sample are 
subsequently recaptured. The summation in the numerator 
tallies the individuals which were in the (kK-r)’ th sample as 
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young (and are thus ‘r’ year old adults at time K) and are 
next seen in some kK or later sample. Thus the expression, 
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of individuals subsequently recaptured, that are ‘r’ years 
old at time K. This proportion of the marked animals, M(k), 
yields the desired estimate. 

This completes the estimation portion of this thesis. 


The remaining estimates are found by substituting into the 


previous equations the estimates for M(k) and N(k,k-r,y,t). 
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